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Density functional theory (DFT) has become the computational method of choice for modeling 
and characterization of carbon dioxide adsorbents, a broad family of materials which at present are 
urgently sought after for environmental applications. The description of polar carbon dioxide (CO2) 
molecules in low-coordinated environments like surfaces and porous materials, however, may be 
challenging for local and semilocal DFT approximations. Here, we present a thorough computational 
study in which the accuracy of DFT methods in describing the interactions of CO2 with model 
alkali-earth-metal (AEM, Ca and Li) decorated carbon structures, namely anthracene (C14H10) 
molecules, is assessed. We find that gas- adsorption energies and equilibrium structures obtained with 
standard (i.e. LDA and GGA) and hybrid (i.e. PBEO and B3LYP) exchange-correlation functionals 
of DFT dramatically differ from results obtained with second-order Mpller-Plesset perturbation 
theory (MP2), an accurate computational quantum chemistry method. The major disagreements 
found can be mostly rationalized in terms of electron correlation errors that lead to wrong charge- 
transfer and electrostatic Coulomb interactions between the molecules. Nevertheless, we show that 
when the concentration of AEM atoms in anthracene is tuned to resemble as closely as possible to 
the electronic structure of AEM-decorated graphene (i.e. an extended two-dimensional material), 
hybrid exchange-correlation DFT and MP2 methods quantitatively provide similar results. 

PACS numbers: 68.43. Be, 73.63.Fg, 81.05.U-, 88.05.Np 



I. INTRODUCTION 

The concentration of carbon dioxide (CO2) in the at- 
mosphere has increased by about 30% in the last 50 
years and is likely to double over the next few decades 
as a consequence of fossil-fuel burning for energy gen- 
eration. This excess of CO2 greenhouse gas may have 
dramatic negative repercussions on Earth's air quality 
and climate evolution. Besides exploitation of renew- 
able energy resources, carbon capture and sequestration 
(CCS) implemented in fossil-fuel energy plants and am- 
bient air have been envisaged as promising cost-effective 
routes to mitigate CO2 emissions [l|. To this end, mem- 
branes and solid sorbents (e.g. activated carbons -AC- 
, hydrotalcites, and coordination polymers -i.e. zeolitic 
imidazolate and metal organic frameworks-) are widely 
considered as the pillars of next-generation CCS tech- 
nologies because of a encouraging compromise between 
large gas-uptake capacities, robust thermodynamic sta- 
bility, fast adsorption-desorption kinetics, and affordable 
production costs Q. A key aspect for the success of 
these materials is to find the optimal chemistry and pore 
topologies to work under specific thermodynamic condi- 
tions. Unfortunately, due to the tremendous variety of 
possible compositions and structures, systematic exper- 
imental searches of this kind generally turn out to be 
cumbersome. In this context, first-principles and classi- 
cal atomistic simulation approaches emerge as invaluable 
tools for effective and economical screening of candidate 
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FIG. 1: MP2 adsorption energy results obtained for a small 
C02/Ca-benzene system and expressed as a function of the 
intermolecular distance z and AO basis set. Ca, C, H, and O 
atoms are represented with large blue, yellow, small blue and 
red spheres, and the dashed line is a guide to the eye. 



carbon adsorbents. 

Density functional theory (DFT) performed with the 
local-density (LDA) and generalized-gradient (GGA) ap- 
proximations of the electronic exchange-correlation en- 
ergy, has become the ab initio method of choice for mod- 
eling and characterization of CCS materials. Standard 
DFT-LDA and DFT-GGA methods are known to de- 
scribe with precision and affordability a wide spectrum 
of interactions in bulk crystals and surfaces, however 
computational studies on the accuracy of DFT in re- 
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producing C02-sorbent forces are surprisingly scarce in 
the literature. In view of the ubiquity of DFT meth- 
ods to CCS science, it is therefore crucial to start filling 
this knowledge gap while putting special emphasis on 
the underlying physics. Computational benchmark stud- 
ies on C02-sorbent interactions, however, are technically 
intricate and conceptually difficult since most CCS ma- 
terials have structural motifs that are large in size. In 
particular, genuinely accurate but computationally very 
intensive quantum chemistry approaches like MP2 and 
CCSD(T) can deal efficiently only with small systems 
composed of up to few tens of atoms, whereas DFT 
can be used for much larger systems (i.e. extended - 
periodically replicated in space- systems composed of 
up to 1,000 — 10,000 atoms). Consequently, computa- 
tional accuracy tests need to be performed in scaled- 
down systems resembling to the structure and compo- 
sition of the material of interest (ejg. the case of organic 
C„H m molecules to graphene [3j-|6|]). Generalization of 
so-reached conclusions to realistic systems however may 
turn out to be fallacious since intrinsic DFT limitations 
(e.g. exchange self-interaction and electron correlation 
errors) can be crucial depending, for instance, on the 
level of quantum confinement imposed by the topology 
of the system. 



In this article, we present the results of a thorough 
computational study performed on a model system com- 
posed of a CO2 molecule and alkali-earth-metal (AEM) 
decorated anthracene (e.g. X-C14H10 with X=Ca, Li), 
that consists of standard DFT (i.e. LDA and GGA), 
hybrid DFT (i.e. PBE0 and B3LYP) and MP2 adsorp- 
tion energy, E a d s 0] , and geometry optimization calcula- 
tions. It is worth noticing that anthracene is structurally 
and chemically analogous to the organic bridging ligands 
found in metal- and covalent-organic frameworks, -MOF 
and COF-, thus our model conforms to a good representa- 
tion of a promising class of CCS materials fsl— fl~2lj . We find 
that both standard and hybrid functionals of DFT dra- 
matically fail at reproducing CO2/AEM-C14H10 interac- 
tions (e.g. E a d s discrepancies of ~ 1 — 2 eV with respect 
to MP2 calculations). This failure is mainly due to elec- 
tron correlation errors that lead to inaccurate electron 
charge transfers and exaggerated electrostatic Coulomb 
interactions between the molecules. In the second part of 
our study we analyse whether our initial conclusions can 
be generalized or not to extended carbon-based materi- 
als, another encouraging family of C0 2 sorbents [TUH. 
For this, we tune the concentration of calcium atoms in 
anthracene so that the partial density of electronic states 
(pDOS) of the model system resembles as closely as pos- 
sible to the pDOS of Ca-decorated graphene. In this case 
we find that DFT and MP2 methods qualitatively pro- 
vide similar results, with hybrid DFT and MP 2 in almost 
quantitative agreement. 
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FIG. 2: Front (F) and top (T) views of equilibrium C0 2 - 
adsorption structures in Ca- and Li-anthracene as obtained 
with standard and hybrid DFT and MP2 methods. Li atoms 
are represented with purple spheres. 



II. COMPUTATIONAL METHODS 



Standard DFT calculations were done using the plane 
wave code VASP [HDl while hybrid DFT and MP2 re- 
sults were obtained with the atomic orbitals (AO) code 
NWCHEM [19]. Numerical consistency between the two 
codes was checked at the DFT-PBE energy level. The 
value of all technical parameters were set in order to 
guarantee convergence of the total energy to less than 
1 meV/atom. Optimized structures were determined by 
imposing an atomic force tolerance of 0.01 eV/A and 
verified as minima on the potential energy surface by 
vibrational frequency analysis. Basis-set superposition 
errors (BSSE) in hybrid DFT and MP2 energy calcula- 
tions were corrected using the counter-poise recipe [20| . 
Indeed, only results obtained in the complete-basis-set 
(CBS) limit HIH can be regarded as totally BSSE free 
however reaching that limit in our calculations turned 
out to be computationally prohibitive due to the size of 
the systems and large number of cases considered. Never- 
theless, we checked in a reduced Ca-benzene system that 
MP2 binding energy results obtained with large Dunning- 
like AO basis sets (i.e. triple zeta cc-pVTZ and quadruple 
zeta cc-pVQZ) and in the CBS limit differed at most by 
20 meV (see Fig.rjJ), thus we assumed the MP2/cc-pVTZ 
method to be accurate enough for present purporses (i.e. 
as it will be shown later, the reported discrepancies are 
in the order of 1 — 2 eV) and regarded it as "gold stan- 
dard". It is worth noting that MP2 results obtained 
with medium and large Poplc-like AO basis sets (i.e. 6- 
311G++(2d,2p) and 6-311G(2df,2pd)) are also in notable 
agreement with MP2/cc-pVTZ results (i.e. E a d s differ- 
ences of 20 — 30 meV in the worse case) whereas MP2/6- 
31G++ estimations (not shown in the figure) turn out to 
be not so accurate. 
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FIG. 3: C02-adsorption energy results obtained for Ca- (a) and Li-anthracene (b), using different optimization and evaluation 
methods. Cases in which optimization and evaluation methods coincide are highlighted with large grey circles. PS, PM, and 
PL notation stands for 6-31G++, 6-311G++(2d,2p), and 6-311G(2df,2pd) Pople AO basis sets, respectively. 



III. RESULTS 

A. AEM-decorated anthracene systems 

Carbon dioxide adsorption energy and equilibrium ge- 
ometry results obtained for AEM- anthracene systems are 
shown in Figs. [2] and [3] Since the number of reaction 
coordinates in C02/AEM-anthracene systems (i.e. in- 
termolecular distances and molecular bond and torsional 
angles) is considerably large, rather than parameterizing 
E a( j s curves as a function of just few of them, first we 
performed DFT and MP2 atomic relaxations and then 
calculated all DFT and MP2 energies in the resulting 
equilibrium structures. In doing this, we disregard local 
minimum conformations, as customarily done in compu- 
tational materials studies, and gain valuable insight into 
the energy landscape provided by each potential. It is 
important to note that due to the size of the systems con- 
sidered we were able to perform tight MP2 atomic relax- 
ations only at the 6-31G++ and 6-311G-r-r-(2d,2p) levels. 
Nevertheless, MP2/6-311G++(2d,2p) and benchmark 
MP2/cc-pVTZ equilibrium structures are very likely to 
be equivalent since results obtained with both methods 
are in fairly good agreement (see Figs. [1] and [3]), and 
MP2/6-31G++ and MP2/6-311G+-(-(2d,2p) equilibrium 
geometries are already very similar (the former case not 
shown here). 



First, we note that equilibrium DFT and MP2 struc- 
tures obtained for both Ca- and Li-anthracene systems 
are perceptibly different (see Fig. [2]). Of particular 
concern is the Ca-anthracene system where, depending 
on the geometry optimization method used, the plane 
containing the CO2 molecule orientates perpendicularly 
(DFT) or parallel (MP2) to anthracene. Also, we ob- 
serve important differences between DFT-PBE and DFT- 
B3LYP optimized geometries, the Ca atom being dis- 
placed towards an outside carbon ring in the hybrid case 
(as it also occurs in the MP2-optimized system). Struc- 
tural differences among Li-anthracene systems are similar 
to those already explained except that the gas molecule 
always binds on top of the Li atom and no off-center 
AEM shift appears in MP2 optimizations. 

Concerning E a d s results (see Fig. [3]), let us concentrate 
first on the Ca-Ci4rli case. As one can see, adsorption 
energies calculated with the same evaluation and geom- 
etry optimization method (highlighted with large grey 
circles in the figure) are very different. In particular, 
DFT methods always predict thermodynamically favor- 
able C02-binding to Ca-anthracene, with DFT-LDA and 
DFT-PBEO providing the largest and smallest E a( j s val- 
ues, whereas MP2/6-3IlG++(2d,2p) calculations show 
the opposite. Moreover, with MP2/6-311G(2df,2pd) 
and MP2/cc-pVTZ methods large and positive adsorp- 
tion energies of ~ 0.6 — 1.0 eV are obtained for DFT- 
optimized geometries, in stark contrast to DFT E a d s re- 
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suits (i.e. ~ —1.0 — 0.0 eV). Adsorption energy dis- 
agreements in Li-Ci4H 10 systems are not as dramatic 
as just described, although the performance of stan- 
dard DFT methods still remains a cause for concern. 
Specifically, DFT-LDA and DFT-PBE predict unfavor- 
able C02-binding whereas hybrid DFT and MP2 meth- 
ods predict the opposite. Also, computed MP2 energies 
in standard DFT-optimized structures are negative and 
noticeably larger than DFT | E a( j s | values. Interestingly, 
the series of binding energies calculated for hybrid DFT 
and MP2/6-311G++(2d,2p) geometries are in remark- 
ably good agreement in spite of the evident structural 
differences involved (see Fig. [2]). 

Since the agreement between MP2 and hybrid DFT 
results in Ca-Ci4Hio is only marginally better than 
achieved with LDA or GGA functionals, common self- 
interaction exchange errors alone cannot be at the root 
of standard DFT failure. Consequently, DFT difficulties 
at fully grasping electron correlations, which in the stud- 
ied complexes account for the 44 % to 57 % of the total 
binding energy, must be the major factor behind the dis- 
crepancy. In fact, upon gas-adsorption important disper- 
sive dipole-dipole and dipole-quadrupole forces appear 
in the systems as a consequence of CO2 inversion sym- 
metry breaking (i.e. the polar molecule is bent), which 
cannot be described by either local, semilocal or hybrid 
DFT approximations. Moreover, the non-linearity of the 
gas molecule also indicates the presence of large electron 
transfers which are well-known to be challenging for DFT 
methods [23[ (i.e. evaluation of excited electronic states 
is involved). But, which of these DFT shortcommings, 
i.e. non-local interactions or charge-transfer errors, is 
the predominant factor behind E a( j s inaccuracies? In or- 
der to get insight into this question, we performed fron- 
tier molecular orbital and charge distribution analysis on 
isolated and joint Ca-Ci4Hi and CO2 complexes [24^ . 
We found that the energy difference between the HOMO 
of Ca-anthracene and the LUMO of the gas molecule, 
A-Efront, varied from 0.7 — 0.9 eV to —1.4 eV when either 
calculated with DFT (standard and hybrid) or MP2/cc- 
pVTZ methods. Positive and large A_Ef ron t values do im- 
ply large reactivity and charge transfers from Ca-Ci4Hio 
to CO2 molecules, AQ. In fact, this is consistent with 
results shown in Fig. [3^, and also with our charge distri- 
bution analysis: AQ values obtained with DFT amount 
to ~ 1 e~, about a 50 % larger than computed with 
MP2/cc-pVTZ. Standard and hybrid DFT approxima- 
tions, therefore, provide overly charged donor (+) and 
acceptor (-) species which in the joint complexes are arti- 
ficially stabilized under the action of exaggerated electro- 
static interactions. In view of this finding, and also of the 
large size of the binding energies reported, we identify the 
inability of DFT methods to correctly describe charge- 
transfer interactions as the principal cause behind its 
failure. Reassuringly, calculations done on Li-anthracene 
systems show a mild improvement on the agreement be- 
tween hybrid DFT and MP2 methods (not seen in the 
standard cases) which is accompanied by smaller A_Ef ront 
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FIG. 4: Partial density of electronic states (pDOS) obtained 
for Ca-anthracene, Ca-graphene, and 3Ca-anthracene using 
the DFT-PBE method. 



discrepancies (i.e. of just few tenths of an eV). Also, 
the amount of charge transferred from Li-Ci4H 10 to CO2 
does not change significantly when either calculated with 
hybrid DFT or MP2 methods (i.e. 0.6 and 0.4 e~, respec- 
tively) . 

Overall, the results presented in this section show the 
key importance of charge-transfer interactions in AEM- 
based CCS nanoparticles which, as we mentioned in the 
Introduction, can be reasonably generalized to similar 
covalent organic structures. Consequently, MP2 or other 
efficient computational schemes embodyin g a lso many- 
electron correlations (see for instance Ref. [25j ) must be 
employed for the rational design and characterization of 
these complexes. 



B. Emulating extended AEM-decorated carbon 
surfaces 



In view of the composition and structure of anthracene, 
it appears tempting to generalize our previous conclu- 
sions to extended carbon-based materials (e.g. AC and 
nanostructures). Nevertheless, we show in Fig. U that 
Ca-anthracene and Ca-graphene are quite distinct sys- 
tems in terms of electronic structure since partially oc- 
cupied s and d orbitals are missing in the former. It 
is important to note that equivalent partial density of 
states (pDOS) dissimilarities are also found when larger 
C„H m molecules are considered (e.g. the case of Ca- 
coronene ||). In consequence of these differences, hy- 
bridization between s<i-sorbent and P-CO2 orbitals lead- 
ing to strong gas attraction will be more limited in Ca- 
C14H10 complexes than in Ca-decorated graphene [13 1. 

Aimed at improving the likeness among the pDOS 
of the model system and Ca-decorated graphene while 
constraining the size of the former system to be small, 
we increased the number of Ca atoms decorating an- 
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FIG. 5: Front (F) and top (T) views of equilibrium (Hem- 
adsorption structures in 3Ca-anthracene as obtained with 
standard and hybrid DFT and MP2 methods. The atomic 
color code is the same than used in previous figures. 
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FIG. 6: C02-adsorption energy results obtained for 3Ca- 
anthracene using different optimization and evaluation meth- 
ods. Cases in which optimization and evaluation methods 
coincide are highlighted with large grey circles. PS, PM, and 
PL notation stands for 6-31G++, 6-311G++(2d,2p), and 6- 
311G(2df,2pd) Pople AO basis sets, respectively. 



thracene [13j, |2fjj. We found that pDOS features in 3Ca- 
C14H10 are already compatible with those relevant for 
CCVbinding to Ca-graphene (see Fig.[4j thus, in spite of 
the obvious chemical and structural deformations intro- 
duced (e.g. now the anthracene molecule is significantly 
bent, see Fig. [5j, we performed additional benchmark 
calculations on this system. 

In Figs. [5] and |6] we show the resulting optimized ge- 
ometries and E a d s values. As can be appreciated stan- 
dard DFT, hybrid DFT and MP2 calculations are now in 
qualitative good agreement: all methods predict equiv- 
alent equilibrium structures and thermodynamically fa- 
vorable C02-binding. Moreover, adsorption energy dif- 
ferences between hybrid DFT and MP2 methods amount 
to less than 0.4 eV in most studied geometries, with hy- 
brid DFT systematically providing the smaller values, 
thus the agreement between these two approaches can 
be regarded as almost quantitative. On the other hand, 



standard DFT approximations tend to significantly over- 
estimate E a( j s (i-e. by about ~ 1.0 — 2.0 eV with respect 
to the MP2/cc-pVTZ method). We note that charge- 
transfer interactions in Ca-overdecorated organic com- 
plexes play also a dominant role and are very intense. 
Indeed, the amount of electron charge transferred from 
3Ca-Ci4rIi to CO2 is AQ = 1.3 e~ according to the 
MP2/cc-pVTZ method and also to hybrid DFT approx- 
imations. In contrast, standard DFT methods predict 
larger AQ values of ^ 1 .6 — 1 .8 e~ . Therefore, we identify 
self-interaction exchange errors as the principal cause for 
overestimation of CO2 binding to carbon-based surfaces 
by standard DFT approaches. 

In the light of the results presented in this section, 
we conclude that standard DFT modeling of extended 
carbon-based CCS materials can be expected to be 
correct only at the qualitative level. On the other 
hand, hybrid DFT approximations conform to a well- 
balanced representation of the relevant interactions in 
AEM-decorated carbon surfaces thus we propose using 
them when pursuing accurate descriptions. 



IV. CONCLUSIONS 

We have performed a thorough computational study 
in which the failure of standard and hybrid DFT meth- 
ods at describing the interactions between X-anthracene 
(X = Li and Ca) and CO2 molecules is demonstrated. 
The origins of this deficiency mainly resides on the in- 
ability of standard and hybrid DFT approximations to 
correctly describe charge-transfer interactions. This find- 
ing has major implications in modeling and characteriza- 
tion of coordination polymer frameworks (e.g. MOF and 
COF) with applications in carbon capture and seques- 
tration. As an effective strategy to get rid of these com- 
putational shortcomings, we propose using MP2 or other 
effective computational approaches incorporating many- 
electron correlations. Moreover, based on the similarities 
in electronic structure found between Ca-graphene and 
3Ca-Ci4rIio systems (and in spite of their obvious chem- 
ical and structural differences), we argue that standard 
DFT modeling of extended carbon-based materials may 
be expected to be correct at the qualitative level. On 
the other hand, hybrid DFT approximations will provide 
quantitative information on these systems. The conclu- 
sions presented in this work suggest revision of an impor- 
tant number of computational studies that are relevant 
to CCS materials engineering. 
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